Measurement-induced qubit state mixing in circuit QED from upconverted dephasing 

noise 
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We observe measurement-induced qubit state mixing in a transmon qubit dispersively coupled to 
a planar readout cavity. Our results indicate that dephasing noise at the qubit-readout detuning 
frequency is upconverted by readout photons to cause spurious qubit state transitions, thus limiting 
the non-demolition character of the readout. Furthermore, we use the qubit transition rate as a tool 
to extract an equivalent flux noise spectral density at / ~ 1 GHz and find agreement with values 
extrapolated from a 1/ f a fit to the measured flux noise spectral density below 1 Hz. 
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High-fidelity measurement is a crucial tool in quantum 
information science. For superconducting qubits [HI2]j 
one widely-used framework for performing quantum non- 
demolition (QND) |3] measurement is the circuit quan- 
tum electrodynamics (cQED) architecture [4j [5] . In 
cQED, a qubit is coupled to a microwave-frequency res- 
onant cavity through a Jaynes-Cummings-type interac- 
tion, in analogy to an atom in an optical Fabry-Perot 
cavity. In the dispersive limit, probing the qubit-state- 
dependent resonant frequency of the cavity implements, 
to first order, a QND measurement of the qubit state. 

In the case of a linear readout cavity [BJ, only re- 
cently has single-shot sensitivity been demonstrated us- 
ing a near-quantum-noise-limited superconducting para- 
metric amplifier [7j |8] , enabling observation of individual 
qubit state transitions in real time [9]. Subsequent ex- 
periments [lUJ [TT] have reported single-shot fidelities of 
94-97%. Nonlinear circuit QED readout methods — using 
either the nonlinearity of the qubit [T^lu"4"| or a nonlinear 
cavity [15] — have shown single-shot fidelities of 86-92%, 
but the former is not QND and the latter is too slow to 
allow continuous qubit monitoring. 

In this letter, we explore the question: How QND is 
a linear cQED readout? We employ single shot read- 
out [H] to directly quantify the rate of measurement- 
induced qubit transitions. We find that dephasing noise 
at the qubit-readout detuning frequency A ro = us q — uj io 
combines with readout photons to induce qubit excita- 
tion and relaxation, making the measurement process no 
longer fully QND. The rate of qubit transitions due to 
such "dressed dephasing" depends linearly on the aver- 
age cavity photon occupation h and the spectral density 
of dephasing noise at the detuning frequency S'(±A ro ), 
consistent with recent calculations which keep higher or- 
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der terms in the dispersive approximation |16j . Further- 
more, the qubit transition rate provides a new probe of 
dephasing noise at |A ro |/27r ~ 1 GHz, a frequency range 
not currently accessible by other techniques. We find 
that our extracted value of dephasing noise at GHz fre- 
quencies is consistent with the "universal" 1/f magnetic 
flux noise [171 I18j typically observed in low frequency 
measurements, suggesting the persistence of this noise 
mechanism over 11 orders of magnitude in frequency. 
Dephasing can be described with the Hamiltonian 

H p = hvf v (t)6- z , (1) 

where f v (t) is a random noise (e.g. flux noise) with 
zero mean and v characterizes the coupling between the 
noise and the qubit. We take vf v (t) to be a small per- 
turbation on the qubit frequency uj q . While the low- 
frequency (<C ujq) components of f v (t) are typically the 
dominant source of qubit dephasing in experiments, the 
frequency spectrum of f v (t) can also have components at 
A ro , which can combine with readout photons to cause 
transitions between the qubit states. The rate for transi- 
tions up and down due to this dressed dephasing is given 

by nana 

r Ti , DD = 4^ 2 S( T A ro )n, (2) 

where g is the qubit-cavity coupling, n is the average cav- 
ity photon occupation, and S(A TO ) is the power spectral 
density of f v (t) at the detuning frequency A ro . This ex- 
pression holds for n <C n cr i t = A 2 D /4g 2 . In the case of 
a symmetric noise spectrum where 5(A ro ) = S{— A ro ), 
we have dd = I^dd- Once the system has reached 
steady state, we expect a spurious excited state popula- 
tion to exist. Using the principle of detailed balance, we 
can express this as 

/ ~ \ i I r t,DD + r t,th _ r t.DD , r t%th 

{a z ) = -1 + « -1 + — h -=r- . (6) 

1 1 + 1 i,DD 1 1 1 1 
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Here l?i = 1/Ti is the intrinsic qubit decay rate, includ- 
ing the Purcell effect, and r Tit h is the qubit's thermal 
excitation rate (the rightmost term here gives the aver- 
age thermal population of the qubit). We have used the 
approximation that T^dd Ii, which is valid for our 
experimental conditions. 

Our experiment, shown schematically in Figure 1, is 
anchored to the mixing chamber of a dilution refrig- 
erator at 50 mK. A transmon qubit [5D] (Ej imax =21.7 
GHz, £,7=220 MHz) is capacitively coupled (g'/2ir = 106 
MHz) to a planar superconducting quasi-lumped-elcment 
readout cavity [HI [22] consisting of a meander induc- 
tor (L=3.7 nH) in parallel with an interdigital capacitor 
(C=175 fF), giving a bare resonant frequency of 6.2724 
GHz. The cavity has asymmetric coupling and is oper- 
ated in transmission; the strongly coupled port sets the 
cavity linewidth k/2tt = 7 MHz. Qubit manipulation and 
readout signals enter from the weakly coupled port via a 
heavily attenuated injection line. Readout photons in the 
cavity acquire a phase shift that depends on the state of 
the qubit, then leave through the strongly coupled port 
and are amplified by a superconducting parametric am- 
plifier (paramp). Four microwave circulators isolate the 
qubit from the strong paramp pump tone. Further am- 
plification is performed by cryogenic and room tempera- 
ture amplifiers (not shown). The output signal is finally 
detected by homodyne mixing and then digitized. This 
method allows us to monitor the qubit state in real time 
and record quantum jumps between qubit energy levels, 
as seen in the inset figure. Further details of the paramp 
mechanism can be found in Refs. [HI H3] & n d supporting 
material. 

In addition, a weakly coupled fast flux line allows mod- 
ulation of the qubit Hamiltonian by noise or coherent 
signals. The fast flux line has a bandwidth of 2.2 GHz, 
defined by a reactive filter at 600 mK and three lossy 
impedance-matched low-pass filters [24] at 600 mK, 100 
mK, and 50 mK. These filters thermalize the line without 
introducing excessive low frequency loss, allowing us to 
pass large currents in the fast flux line without heating 
the mixing chamber. To calibrate the coupling of the fast 
flux line to the qubit loop, we extract the flux-to-qubit- 
frequency transfer function from qubit spectroscopy. We 
then measure the qubit frequency as a function of ap- 
plied dc current through the fast flux line using Ram- 
sey fringes. Combining these factors with the measured 
frequency-dependent attenuation of the fast flux line al- 
lows us to convert room-temperature power into a flux 
in the qubit loop. The coupling is sufficiently weak (120 
mA/$ ) that Johnson noise from the 50 fl impedance of 
the fast flux line is not the dominant source of flux noise 
in the qubit loop for frequencies at or below w g /27r. 

For each qubit bias point, the dispersive shift 2\ is cal- 
culated using an expression that accounts for the higher 
excited states of the transmon qubit and includes correc- 
tions for the qubit-induced Kerr nonlinearity |13j . We use 
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FIG. 1. Experimental setup and readout trace (color online). 
The transmon qubit is coupled to an asymmetric microwave 
readout cavity. Noise or coherent tones can be injected into 
the qubit loop via a weakly coupled fast flux line with 2.2 GHz 
bandwidth (grey dashed lines indicate flux coupling). The 
readout signal is amplified by a superconducting paramp, en- 
abling continuous high-fidelity monitoring of the qubit state. 
The inset shows a sample output trace (red line) and the cor- 
responding extracted qubit state (dotted black line). 



the dispersive shift information to calibrate the number 
of photons in the readout cavity using the ac Stark shift 
[25] and to choose the readout frequency ui IO /2ir. We se- 
lect a readout frequency halfway between the cavity reso- 
nant frequencies corresponding to the qubit in the ground 
and first excited states u IO — ^ [w cav (|0)) + 
w C av(|0)) + X- This choice of readout frequency means 
that the average cavity photon occupation n is unaffected 
(up to a 5-10% correction at the highest n used in the 
experiment) by whether the qubit is in state |0) or |1), 
simplifying the analysis of the dependence of transition 
rates on n. Qubit coherence times varied monotonically 
depending on the qubit frequency, from T\ — 290 ns and 
T 2 = 550 ns at w 9 /2vr = 5.705 GHz to T x = 910 ns and 
T 2 = 1.35 /is at ui q /2ir — 5.075 GHz. These numbers 
represent T\ values about a factor of two below the Pur- 
cell limit, and pure dephasing times T v much longer than 
Ti. 

The measurement protocol consisted of readout pulses 
lasting 17.5 fis occurring every 100 fis. The long de- 
lay ensured that the qubit would fully relax to its ther- 
mal ground state (~1.4% excited state population, cor- 
responding to a qubit temperature of 60 mK) between 
measurement runs. We took 10 4 individual time traces 
for each combination of experimental parameters. The 
measurement traces were analyzed by smoothing to op- 
timize signal-to-noise ratio and then using a hysteretic 
threshholding algorithm similar to that demonstrated in 
[2"B] to determine the qubit state at each time point. This 
allowed qubit state populations to be determined both at 
steady state and as a function of time into the readout. 
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FIG. 2. Spurious excitation with coherent fast flux tone (color 
online). The qubit steady state population during measure- 
ment with a coherent microwave tone applied to the fast flux 
line is shown ELS cLS cl function of ujg and n. The panels cor- 
respond to RMS added fluxes of 300 fi&o (a) and 550 fi&o 
(b). The black x's denote the values of A ro (n) extracted in- 
dependently from qubit spectroscopy. Regions with no data 
(at the lowest and highest frequencies) are shown in uniform 
dark blue. Panel (c) shows the qubit population versus ujs 
for n = 6 and an RMS added flux of 420 fi&o (black dotted 
line and diamonds) , along with the corresponding qubit spec- 
troscopy trace (red solid line). The horizontal axes are aligned 
such that the detuning between readout and spectroscopy fre- 
quencies is the same as the fast flux frequency. Panel (d) 
shows the center (solid line) and widths (dotted lines) of the 
qubit excited state population response for the data in (b). 
The red (bold) traces are experimental data, while the blue 
(narrow) traces are from numerical simulations. 

The population extraction algorithm was tested on simu- 
lated data traces with realistic experimental parameters, 
and shown to give fractional errors of less than 5% for 
population estimates [23"] . 

To test the dressed dephasing theory, we began by in- 
jecting a continuous microwave tone into the fast flux 
line at a frequency tug ~ |A ro (n)|, where we note that 
A ro depends on n due to the ac Stark shift [35]. This 
tone produces a small flux oscillation in the qubit loop; 
the power Pg of the tone was varied to produce RMS 
flux excitations of up to 825 /i'&oi corresponding to qubit 
frequency fluctuations of up to 3-6 MHz RMS (depend- 
ing on qubit bias point). These fluctuations are much 
smaller than the ac-Stark-broadened qubit linewidth for 
the values of n studied. For each value of n, we stepped 
the ug through a range of around 100 MHz centered on 
|A ro (fi)|. 

The results of this experiment (with the qubit biased 
at Lo q /2ir — 5.075 GHz) are shown in Figure [2] When 
Pg ^ 0, qubit state mixing occurs as long as u>g is within 
roughly a qubit linewidth of the detuning frequency, and 



is most noticeable when u)g =|A ro (n)|. Panels (a) and 
(b) show the qubit excited state population with vary- 
ing n for P ff =300 /i<I>o and 550 n&o, respectively. The 
value of |A ro (n)|, found independently from qubit spec- 
troscopy, is denoted with black crosses. Figure [2jc) plots 
the spurious excitation as a function of ujg along with 
the independently measured spectroscopy signal for the 
same n; the horizontal axes are aligned such that the 
fast flux frequency is the same as the detuning between 
the readout and spectroscopy frequencies. This choice 
highlights the correlation between the values of ujg which 
cause qubit state mixing and A ro . 

Figure [2^d) shows the extracted center frequency and 
linewidth of the regions of qubit excitation in Fig. [2^b). 
We also plot the center frequency and linewidth of qubit 
excitation, for the same experimental parameters, cal- 
culated from numerical simulations of the multilevel 
Jaynes-Cummings Hamiltonian with an added coherent 
flux tone. The simulations agree well with the experi- 
mental data on the location and width of the peak, with 
some width discrepancy appearing at high n. The simu- 
lations also agree to within a factor of ~ 1.4 or better (de- 
pending on bias point) on the steady-state qubit excited 
population. We attribute the remaining discrepancy to 
uncertainty in the calibration of our fast flux line. 

In the absence of a fast flux tone, we still observe 
some spurious qubit excitation with increasing n, which 
we postulate is due to the intrinsic qubit flux noise at 
A ro (n) being upconverted by readout photons. The ef- 
fect is about 1% additional excited population at steady 
state per 10 photons cavity occupation, increasing more 
rapidly at higher photon numbers (up to 10-15% excited 
state population for n m 40) . This measurement-induced 
state mixing, which reduces the fidelity of cQED mea- 
surement, can also be used as a spectrometer for dephas- 
ing noise at A ro . 

We examined this notion by intentionally applying flux 
noise to the qubit loop using the fast flux line and observ- 
ing spurious excitation during measurement. The noise 
was generated by amplifying the Johnson noise of a room- 
temperature 50 VL termination with variable gain. The 
experiment was performed with white noise filtered to lie 
either in the band from 10 MHz to 2.2 GHz or from 180 
MHz to 2.2 GHz. The steady-state qubit populations 
were essentially identical between these two types of ap- 
plied flux noise, suggesting again that only dephasing 
noise components near A ro (ri) are responsible for spuri- 
ous excitation. 

Figure |3][a) shows the steady-state qubit excited popu- 
lation as a function of n and the spectral density of added 

1 /2 

flux noise at the detuning frequency add (A ro (n)). The 
qubit excited state population scales roughly linearly in 

1/2 

n and quadratically in S' $ / add (A ro (n)), as predicted by 
the dressed dephasing theory. To ensure that the added 
noise is not causing excitation in the absence of measure- 
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FIG. 3. Qubit population with added noise (color online). 
Part (a) shows the steady-state qubit population as a func- 
tion of added flux noise at |A ro (n)| and the number of mea- 
surement photons n in the readout cavity. Part (b) shows 
the qubit population immediately after the readout has fully 
energized. A thermal population of 1.4% is visible, and is 
equivalent for all values of added flux noise and n. 



ment, we also examine the initial qubit population at the 
start of measurement from the same data set, shown in 
Figure^b). We find that the excited state population at 
the start of the measurement corresponds to the thermal 

1/2 

population, and is independent of n and add (A ro (n)). 

Given the very weak charge dispersion of the trans- 
mon [20 and the low level of critical current noise in 
similar sub-micron Josephson junctions [27] . we expect 
flux noise to be the dominant source of dephasing in our 
qubit sample [25]. If we attribute all spurious excita- 
tion to upconverted flux noise, we can extrapolate the 
data in Fig. [3ja) back to a value of the added flux noise 
where there is no excitation above the thermal popula- 
tion, yielding an estimate of the intrinsic flux noise at 
the detuning frequency. Using this method, we extract a 
flux noise spectral density 5^ /2 (705 MHz) = 0.009±0.004 
/i^o/VHz. Using experimentally measured values for g, 
A ro , v, T\, and the slope of (<j z ) versus n, the theoretical 
dressed dephasing expressions Q and ^ give a value of 
Sl /2 (705 MHz) = 0.011 ± 0.002 /x$ /VHz. 

We performed two other experiments to measure flux 
noise at other frequencies. First, we measured Ramsey 
fringes at a rate of one per second for 90 minutes, then 
fit each fringe to extract the qubit frequency. The fluctu- 
ations in the qubit frequency over time can be translated 
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FIG. 4. Spectral density of flux noise vs frequency (color on- 
line). We plot equivalent flux noise extracted from Ramsey 
fringes (black line), Rabi oscillation decay (red circles), and 
our measurements of spurious excitation (blue squares). The 
red line is a fit to the Ramsey fringe data using a 1 / f a power 
law, with the grey shaded area representing the 95% confi- 
dence interval on the fitted value of a — 0.57. The inset view 
shows a detail of the Ramsey fringe data and fit line. 



into an effective flux noise spectral density for frequen- 
cies below 0.5 Hz. Second, we measured Rabi oscillations 
at different Rabi frequencies £Ir and extracted the decay 
rate r 2 . We used this decay rate to obtain a spectral 
density of qubit frequency fluctuations Sg^ at FIr using 

the relation T 2 = fTi + ~T U , where T v = irSsu q (^R) 
[2"9"] . Converting Sg u (£Ir) to an effective flux noise gives 
us data for frequencies between 1 and 20 MHz. 

The flux noise values extracted by these methods, as 
well as those from our spurious excitation data, are shown 
in Figure [4] We fit the Ramsey data to a 1// Q power 
law [T7J [3U] corresponding to the red trend line. This 
fit agrees with the extracted values from both Rabi de- 
cay and spurious excitation data, representing a power 
law for flux noise that appears to hold over 11 orders 
of magnitude in frequency. The fit coefficients give 

-rl/2 



a = 0.57 ± 0.03 and S% (1Hz) = 1.5 ± 0.1/i$o/VHz, 
both of which agree with typical values reported in the 
literature [TTJ [3TT435J . We note that other recent work 

1 /2 

produced a similar (1Hz) but found a = 0.9 — 1.0 
[301 136j ; sample-to-sample variation in a of this magni- 
tude has been noted elsewhere [TTJ [M] [35] . 

The correspondence of the low-frequency fit to the ex- 
tracted flux noise at A ro (n) suggests that the "universal" 
low-frequency flux noise |18j persists to GHz frequen- 
cies. In addition to dephasing the qubit state, this noise 
also reduces the QND character of circuit QED measure- 
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ment, in agreement with the dressed dephasing theory. 
This sets limits on achievable readout fidelity — even with 
quantum-limited post-amplification — by giving a penalty 
for increasing n. Recent efforts to understand and im- 
prove low-frequency dephasing noise [551 133 EE] m ay 
therefore also provide a route to improved qubit read- 
out fidelity. 
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